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This paper presents an analytical method to determine static and dynamic char- 
acteristics of annular parallel-grooved seals. The governing equations were 
derived by using the turbulent lubrication theory based on the law of fluid 
friction. Linear zero- and first-order perturbation equations of the govern- 
ing equations were developed, and these equations were analytically Investi- 
gated to obtain the reaction force of the seals. An analysis Is presented that 
calculates the leakage flow rate, the torque loss, and the rotordynamlc coeffi- 
cients for parallel-grooved seals. To demonstrate this analysis, we show the 
effect of changing number of stages, land and groove width, and Inlet swirl on 
stability of the boiler feed water pump seals. Generally, as the number of 
stages Increased or the grooves became wider, the leakage flow rate and rotor- 
dynamic coefficients decreased and the torque loss Increased. 


INTRODUCTION 

Annular pressure seals can significantly Influence the dynamic behavior of 
rotating machinery by the presence of a high-pressure difference In the close 
clearance spaces of the leakage path. 

Black and Jenssen (refs. 1 to 3) have explained the Influence of seal 
forces on the rotordynamlc behavior of pumps. Childs has analyzed the short 
seal (ref. 4) and has made finite-length analyses (ref. 5) based on Hlrs 1 gov- 
erning equation, which yields an analytical expression for the seal dynamic 
coefficients Incorporating all of Black and Jenssen's various developments. 

Although these results apply only for small seal motion about a centered 
position, Allaire, et al. (ref. 6), the authors (ref. 7) have expanded these 
analyses to calculate dynamic coefficients at large eccentricities. 

Fleming (ref. 8) has developed an analysis for gas seals with a constant 
clearance or with convergently tapered geometries. Child (ref. 9) Investigated 
dynamic coefficients for convergently tapered seals both analytically and 
experimentally. 
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Previous analytical and experimental developments have generally examined 
dynamic characteristics of the annular straight seal and the tapered seal. But 
the dynamic characteristics of the parallel-grooved seal have not been analyzed 
theoretically. Bolleter (ref. 10) experimentally Investigated stability limits 
for balance pistons with two different types of serration. He showed that ser- 
rations that are deep and wide prevent seizure. 

In this paper the static and dynamic characteristics of the multistage 
parallel-grooved seal operating within the turbulent flow region are analyzed 
theoretically with consideration of the Inertia effect. The present analysis 
combines the previous analysis of the straight seal with the analysis of the 
labyrinth seal performed by the authors. Namely, land analysis Is used for the 
straight seal and groove analysis Is used for the labyrinth seal. 


SYMBOLS 
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nominal seal radial clearance, cm 

seal damping coefficients, N s/m 

journal diameter, cm 

fluid force, N 

seal radial clearance, cm 

seal stiffness coefficients, N/m 

seal length, cm 

land width, cm 

groove width, cm 

2 

seal add mass coefficients, N s /m 
fluid mean depth, cm 
fluid pressure, MPa 

3 

rate of leakage, m /s 

seal radius, cm 

axial Reynolds number 

circumferential Reynolds number 

number of stages 

groove depth, cm 

time, s 
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Torq torque, N m 

u,w tangential and axial fluid velocity components, m/s 

V journal surface velocity, m/s 

x circumferential coordinate 

y radial coordinate 

z axial coordinate 

c small eccentricity ratio 

0 diverging angle of stream behind land section 

x friction loss coefficient 

fluid viscosity, mPa s 

2 

fluid kinematic viscosity, m /s 
\ loss coefficient 

3 

P fluid density, kg/m 

t shear stress, Pa 

« journal angular velocity, rad/s 

Subscripts: 

a axial direction 

c radial direction 

d groove defined In equation (13) 

ex exit 

f between clearance flow and cavity flow 

g groove 

In Inlet 

j journal 

1 land 

r circumferential direction 

s casing 
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GOVERNING EQUATION AND ANALYSIS 
Governing Equation 

Figure 1 Illustrates the geometry of the parallel-grooved seal. Under the 
usual assumptions for problems of throughflow across annuli with a fine clear- 
ance, the momentum and continuity equations are, respectively, as follows: 
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The fluid velocities u and w are Integrated across the film and are 
transformed Into the mean velocities u m and w m . 
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For a fully developed turbulent flow regime the velocity profile shape becomes 
flat, and the quantity Is close to unity as shown In Burton's experiments 
(ref. 11). 
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Equations (1) to (3) can be rewritten as 
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Assumptions for Analysis 

The assumptions for this analysis are as follows: 

(1) The fluid Is liquid and Incompressible. 

(2) The fluid flows Into the groove chamber with a constantly diverging 
angle (ref. 12). This Is Illustrated In figure 2. 

(3) For a small motion of the seal journal about a centered position, the 
streamline In the groove moves with the journal. 

(4) The groove cross section Is rectangular. 


DERIVATION OF SHEAR STRESS IN MOMENTUM EQUATION 

The shear stress terms of equations (4) and (5) are discussed here. 

Land . - The axial components T £sa> T zja of shear stress at the 
casing and journal are given by 

Isa T l j a T la 2 ^^la 3 Im 
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where \^ a defines the friction coefficient between the flow and the wall 
surface: 
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The circumferential components T£sr* T £jr of shear stress at the 
casing and journal are given by 
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( 10 ) 
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This equation Is strictly applicable only for R a » R r . Therefore the shear 
stress of the land can be written as follows: 
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Where the effective viscosity y a (ref. 13) Is 

u » 0.01(7- 3&)R A. M 0 

a ala 



Groove . - The crossflow In the groove Is considered to be the cavity flow 
plus the clearance flow (fig. 2). The exchange of energy within a small mixing 
area between both flows Is Influenced by the entrance velocity and the geo- 
metric shape of the groove. The clearance flow can be described by a stream 
tube If this mixing area Is replaced by a separating layer. The cavity flow 
assumes that the momentum that Is supplied from the journal Is balanced by the 
sum of the momentum lost by fluid friction at the separating layer and the 
momentum due to cavity flow. If the cavity flow Is replaced with flow through 
a circular pipe, the fluid mean depth of groove MG Is given as 
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The friction coefficient between the groove wall and the fluid Is 
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Therefore the momentum that Is supplied from the side and root of the qroove 
Is expressed, where 


i px d 


( R - Ti - 4J rh- u dm } 2 * ( T+(T-L g tan6 ) +lJ 


(14) 
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And the momentum that Is lost from the fluid friction between the clear- 
ance flow and the cavity flow is 


-u ) 2 - L /a osd 
2 f ' dm gm 1 g 


(15) 


where \f defines the fluid friction loss coefficient between the cavity flow 
and clearance flow, and \f = 0.1 (ref. 14). Using equations (14) and (15), 
one can calculate the circumferential fluid velocity of the cavity flow Uq m . 
The axial component of shear stress Is given by 
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Therefore the shear stress of the groove can be written as 
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Derivation of Static Characteristics 

Axial fl u id velocity. - The pressure loss is stated as follows: 
The inlet loss at the land entrance is 


A P lin ~ 2 P ^l U lo 
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(after 2nd stage) 
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The wall friction loss In the land Is 


LP l ■ * -^ T "ha U U 


The exit loss due to diverging flow behind the land Is 


LP lex ”Y^2 U lo 


where t* 16 exit loss coefficient 


" ( J ~c, +r~taZQ ) (at each stage) 
' lo g ' 

^2 = [ 1 ~cT L Tt'] (at seal exit) 
lo 


The friction loss at the wall and the separating layer of the groove is 

/ C \ 2 ^ 

The pressure restoration due to deceleration at the groove Is 
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The pressure drop across the seal Is equal to the sum of each pressure 
loss and Is stated as follows: 
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Rate of leakage . - The main advantage of the grooved seal Is that leakage 
flow can be minimized but seal components need never rub. Seal leakage flow 
may also be denoted by Q, where 

fR+T 
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Figure 3 shows the calculated results of leakage flow. As the number of 
stages Increased or the land became narrower, the leakage flow decreased. 

Circumferential fluid velocity . - From the momentum equation (5) for the 

land 
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where z = C ^ 0 /x - a and V = Rw. The boundary condition Is z = z n _-|/ 2 ; 
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The circumferential fluid velocity u^ Q at the land Is 
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From the momentum equation (5) for the groove 
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where u<j m can be calculated from equations (14) and (15) and the boundary 
condition Is 


a - *. 


V = U 9o(z n > 


The circumferential fluid velocity u go at the groove Is 
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Figure 4 shows the axial distribution of circumferential fluid velocity 
AP = 0.49 MPa, N = 4000 rpm, L = 55 mm, and S = 20. The results Indicate that 
the circumferential velocity approaches one-half of the shaft angular velocity 
exponentially. 

Torque loss . - The method for estimating the torque loss In a grooved seal 
treats the loss as frictional dissipation by viscous shear In an annulus. From 
equation (10) the torque loss of the land Is 
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The torque loss of the groove assumes that the groove Is divided Into four 
parts as shown In figure 5. The shear stress of part 1 Is 
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Therefore the torque loss of part 1 Is given as follows: 
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Therefore the torque loss of the groove Is 
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And the total torque loss of the parallel-grooved seal is defined by the fol- 
lowing equation: 


S+l s 

AT = Z2 Torn 1 + XL Torq (33) 

t=2 L i=l 9 


Figure 6 shows the torque loss of the parallel-grooved seal for 
AP = 0.49 MPa and L = 55 mm without inlet swirl. As the number of stages 
increased or the seal clearance decreased, the torque loss increased. 


Derivation of Pressure Distribution 

For a small motion about a centered position the clearance, pressure, and 
velocity are expanded in the perturbation variables as follows: 

u m C + e<|> P = P + eP, V - w + eu, »<_•«* eu ( 34 ) 

“ o o 1 mo 1 mo 1 

where C 0 , P 0 , w 0 , u 0 are steady-state values and y, P-| , w-| , u-| are small 
perturbations. The short-bearing solution (ref. 4) is developed for the first- 
order equation by neglecting u-| , the pressure-induced circumferential veloc- 
ity component. 

Static pressure distribution . - The steady-state equation described is 
static, has zero eccentricity, and is solved analytically. Substitution of 
equation (6) into equations (11), (18), and (34) yields steady-state, axial- 
direction momentum equations for the land 
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and for the groove 
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where \ sa defines the equivalent friction coefficient of the groove. 

A = A +0 . 2SX .~2tanQ 
a a ga f 

Pressure distribution of nonsteady state . - Substitution of the perturba- 
tion variables of equations (6) and (7) yields the following perturbation equa- 
tion for the land: 
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Equation (38) can be Integrated by using the boundary condition z = Zn-1/2* 
w Zl = w ZKz n _ 1/2 ) as follows: 
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Equation (19) yields the following perturbation-variable boundary condi- 
tion at the seal Inlet: 
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The perturbation-variable boundary condition at the seal exit Is given as 
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The complete solution for the perturbation pressure Is obtained as 
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where PRL1 to PRL8 are provided In appendix A. 


Similarly for the groove the axial -momentum equation and the continuity 
equation reduce to 
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Substitution of equation (29) and the boundary condition z = z n , w g -| = w gl ( Z ) 
Into equation (44) yields n 
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The pressure boundary condition Is 
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The perturbation-variable boundary condition Is obtained as 


2-2 


(46) 


111 


(47) 


gl(z n ) 2 W lo W gl( 2^7 


The complete solution for the perturbation pressure Is obtained as 
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where PRG1 to PRG2 are provided In appendix B. 


Dynamic Force 

For a small motion about a centered position the clearance function Is 
defined In terms of the radial seal displacement (AX, AY). 


Eip = - AT aos$ - AY sin$ 
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The components of the reaction force acting on the seal journal are defined by 
the Integrals 
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Substituting for and Its derivatives Into equations (42) and (48) yields 
the following form for the seal coefficients: 
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Because the seal coefficients of each stage are different, the following 
definitions are used. The dynamic coefficients of the parallel-grooved seal 
become 
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The dynamic seal coefficients for the land are 
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where IPRL1 to IPRL8 are provided In appendix C. The dynamic seal coefficients 
for the groove are 
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where IPRG1 to IPRG21 are provided In appendix C. 


NUMERICAL EXAMPLE 

Seal coefficients were calculated for a pump seal with the characteristics 
shown In table I. The dynamic coefficients decreased with an Increase In the 
number of stages (fig. 7); that Is, the axial Reynolds number became small 
because the pressure loss Increased with an Increase In the number of stages. 
The cross-coupling terms were relatively sensitive to swirl at the seal 
entrance; as the Inlet swirl was propagated right through, the mean circumfer- 
ential fluid velocity approaches one-half the journal speed exponentially 
(ref. 15). 

Figure 8 Illustrates the Influence of the ratio L 7g of land width L ? 
to groove width plus land width (L, + L g ) for N = 4000 rpm and L/D = 0.25. 

As the clearance ratio C^/Cg became small, the dynamic coefficients became 
large with an Increase In^the ratio L^g. 


CONCLUSIONS 

The static and dynamic characteristics of the annular parallel-grooved 
seal were theoretically Investigated with consideration for the effect of the 
turbulent flow and the inertia term. 

As the number of stages Increases or the land becomes narrower, the 
dynamic coefficients and the leakage flow rate decrease. But torque loss 
Increases rapidly. 
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TABLE I. - NUMERICAL CALCULATION MODEL 


Working fluid Water 

Fluid temperature, K 293.15 

Density, P , kg/m3 9.982x102 

Viscosity, v, mPa s 1.009 

Kinematic viscosity, v, m2/s 1.006x10-6 

Journal radius, R, mm 100.0 

Seal radial clearance, C 0 , mm 0.4 

Seal length, L. mm 55 to 205 

Groove depth, T, mm 3.0 

Divergent flow angle, e, deg 4.0 

Journal rotating frequency, N, rpm 2000 to 8000 

Pressure difference, AP, MPa 0.49 to 4.9 
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Figure 1. - Geometry of parallel-grooved seal. 



Figure 2. - Streamlines and coordinate system for seal analysis. 
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CIRCUMFERENTIAL VELOCITY Uo (n»/s) 



Figure 3. - Leakage flow rate. 



Seal Profile 


Figure 4. - Axial distribution of circumferential fluid velocity 
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